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Abstract 

We study the role played by the final-state interactions (fsi) in heavy quark 
decays using Pomeron and Regge exchanges to describe high energy scatter- 
ing. At center of mass energy \/s ~ 5GeV , Pomeron dominance does not 
apply. We study the behavior of the decay amplitudes as \Js is varied close 
to the B mass. We also investigate the behavior of the decay amplitude as 
y/s — * oo. Our conclusion is that the decay amplitudes approach a real value 
asymptotically. 
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Recently, much has been written |l]-[5] about the role played by the final state interactions 
(fsi) in heavy quark decays. The conclusions are contradictory: Ref || claims that fsi 
continue to play an important role in B decays by imparting non-vanishing phases to the 
decay amplitudes ( a conclusion the authors of |||| exploit ) while Refs [|1| and || arrive at 
an opposite view that the phases generated by fsi are insignificantly small. 

In the present work, we have analyzed the question anew to achieve a better understand- 
ing of the role of fsi at B mass and higher. We begin with a few relevant details which also 
assist in establishing the notation. The first point to realize is that at B mass, center of 
mass energy yfs — tub, the total cross-sections for the most hadronic reactions (K + p, K + n 
excepted ) are decreasing with energy ||. Thus Pomeron dominance does not set in at 
y/s ~ 5GeV. Regge exchanges play an important role at B mass, though there is an energy 
at which Pomeron dominance sets in signaling the rise in total cross-sections. This occurs 
at energies above ~ 15GeV. 

In order to follow the reasoning of , define the partial- wave decomposition of scattering 
amplitude, T(s, t), for dissimilar spin-less particles as follows, 

T{s,t) = ^J2( 2l + 1 )Ms)Pi(cos9), (1) 

where k is the center of mass momentum and Ai(s) the partial- wave amplitudes. The 
discontinuity for the two-body decay amplitude, M(s), for B —> f, where f involves spin-less 
particles, arising solely from the elastic channel is 

a(s) = DiscM(s) = i^L_ / j^^HV-s - * - E,yr{., t)M(s), (2) 

where T(s,t) is the amplitude for / — > f scattering. Using (1) in (2) one gets 

a(s) = A* (s)M(s). (3) 

In what follows we ignore the spin and isospin complications as they are inessential to our 
arguments. Because of this only A (s) will appear in our formulation. If M(s) is a real 
analytic function satisfying a dispersion relation, the spectral function a(s) is real. For a 
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single-channel problem, the S-wave amplitude A (s) is e iS °sin5o. Thus the decay amplitude 
in (3) must have a phase 5q in order to generate a real cr(s). This is the content of Watson's 
theorem 0. In general, when other channels are present, Aq(s) = {rje 2lSl — l)/2i, where Si 
is an eigenphase and the elasticity is limited by < rj < 1. Reality of a(s) now requires 
additional contributions to the right hand side of (3) which we discuss in the following. 

To be specific we consider the decay B — > D*ir but ignore the spin and isospin of the final 
state particles only to the extent needed. As noted in f3j, Dp channel will mix with D*ir, but 
D*7r' and Dp' channels, where tt' and p 1 are higher mass particles with the same quantum 
numbers as tt and p respectively, could also mix with D*n state. Indeed, intermediate states 
D*(nir) enter the picture, where for n odd one could form a 7r-like and for n even a p-like 
combination. A generalization of (3) for the spectral function is 

a(s) = DiscM B -,D»M = A* (D*n -> D*ir)M B -* D *M + A o( D *P ^ D*ir)M B -*D* P (s) + 

A*(D*n' -> D*ix)M B ^ D ^{s) + A* (Dp' -> D*ir)M B -> D A s ) + ( 4 ) 

Let us now assume that a Regge-exchange description applies to the scattering ampli- 
tudes. (Ref lU has concluded that cuts do not play an important role. ) The elastic channel 
D*ir — > D*tt receives contributions from the Pomeron (P), and the p— and /— Regge tra- 
jectories. The inelastic scatterings Dp — > D*tt and Dp' — > D*n receive contributions only 
from 7r or tt' Regge trajectories. However, the inelastic scattering D*ir' — > D*tt, as well as 
Dp — > Dp', receives contributions from P, f, p and p' trajectories. Let us now study the 
S-wave amplitudes and their behavior ClS db function of s. 

The Pomeron contribution to the elastic channel, D*n — > D*tt , or the inelastic channel 
D*7r' — > D*tt, can be parameterized in a form (we assume that Pomeron amplitude is 
absorptive) 

T p (8,t)=ij3{±) 1M J"*, (5) 
So 

where the Pomeron trajectory is represented by ap(t) = 1.08+a' P t, bp = a' P ln(-^), a' P = 
0.25GeV~ 2 , and t m — s(l — cos6)/2 ( s and t are expressed in GeV 2 with s = lGeV 2 ). 
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Projecting the S-wave amplitude from (5) in the limit where all masses can be ignored, one 
obtains 

IbubpSQ So 

Typically, for (3 « 10 A£(s) « 0.3* at ^/s = SGeF. The important point is that A%(s) 
is purely imaginary and small compared to unity Because of the s dependence of bp in 
(6), Aq(s) decreases slowly with s along the imaginary axis of the Argand diagram, Fig.l, 
reaching a minimum, of approximately half of its value at y/s = 5GeV, at \fs ~ 500GeF. 
Aq(s) then increases slowly eventually violating the unitary limit, ImAi(s) < 1, at y/s ~ 
2.5 x 10 u GeV for (3 = 10. This, of course, implies that op(0) cannot remain at 1.08 for all 
energies, and must eventually decrease. Ref. || makes the argument for "softening" of the 
Pomeron due to the violation of the Froissart bound on the total cross-section albeit at 
an energy several orders of magnitude higher. Around the B mass, say, \fs ~ 5 to lOGeV^, 
the Pomeron partial-wave amplitude, Aq, is purely imaginary and decreases with energy. 

Consider now the Regge-exchange amplitudes. As a typical example, the p-exchange 
amplitude [|[] is represented as, 

1 {acp(t)) sirnia p [t) s 

For simplicity, assume a p {t) = 0.5 + a't, with a' ~ \GeV~ 2 and f| T(a p {t))sin(iTa p (tj) « 
r(a p (0))sm(7ra p (0)) ~ y^. We then obtain the S-wave amplitude as 

< P) ( S ) « o/ff, (-r (0) (0-30 + ^0.32). (8) 

For /3(0) ~ 100 |Q and at -^/s ~ 5GeV, one finds, first, that this amplitude is much smaller 
than the Pomeron-generated amplitude, the real and imaginary parts being « 0.035 and, 
second, that A p \s) decreases with energy as (~) Qp< ' ' ) Thus in the Argand diagram 

of Fig.l, while the Pomeron amplitude moves slowly downwards along the imaginary axis 
with increasing energy around the B mass, the Regge amplitudes move towards the base 
("no scattering" limit ) of the unitary circle much faster. 



Im A, 



Re A, 

-0.5 0.5 

Fig.l The partial- wave amplitudes in Argand diagram. P represents the Pomeron-generated, and R 
the Regge-generated amplitude. The arrows indicate the motion of the amplitudes with increasing 
s (for ^/i < 500GeV). 

The foregoing analysis allows us to understand the conclusion reached in Refs [|]] and 
||. Define the I = S-matrix as ( P and R stand for Pomeron and Regge respectively), 

S I=0 (s) = l + 2i(A£(s) + A*(s)). (9) 

As discussed above, ImA^s) and ImA^s) are both small but positive. Thus, in the B 
mass region, the real parts of the diagonal elements of S are close to, but less than, unity. 
The imaginary parts of the diagonal elements, arising entirely from Aq(s), are necessarily 
small. In the schemes of [[jj and [Q], the off-diagonal elements of the S-matrix arise entirely 
from Aq(s) and, therefore, are small. As a consequence, the two-channel unitarity is not 
satisfied. A difference between our formulation and those of QTJ] and |3j is that we include 
off-diagonal channels, D*tt' — > D*tt and Dp' — > Dp etc., which are Pomeron driven and 
could, in principle, dominate the off-diagonal elements of the 5*-matrix. 

As we increase the energy but still keeping it in the B—mass region, the real parts of the 
diagonal elements of the S-matrix rise towards unity and the imaginary parts decrease. The 
coupled channel problem involving channels D*n, D*n' , D*tt", Dp, Dp', Dp", tends 
towards a block- diagonal form with D*n, D*n' , D*n" states forming one block and 
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Dp, Dp' , Dp", composing the other. The transitions between these blocks, which proceed 
only through Regge exchanges, get weaker with increasing energy. The transitions within 
each block being Pomeron dominated result in largely real S-matrix elements (if the Pomeron 
is purely absorptive). 

The question: what happens in the limit s — > oo? is harder to answer. As we saw, this 
limit cannot be taken for the Pomeron amplitude with «p(0) = 1.08. However, if we go to 
an energy large yet smaller than the value where violation of partial-wave unitary sets in, 
the Regge-driven Aq(s) approach the base of the unitary circle of Fig.l and the Pomeron- 
driven amplitudes remain on the imaginary axis. The resultant S-matrix approaches a 
real-symmetric form. If, further, the Pomeron "softens" to the extend that Op(0) — > 1 
asymptotically, then Aq (s) also approaches the base of the unitary circle and the S-matrix 
approaches the unit matrix ("no-scattering" limit). 

Let us return to the case of B decays and study the energy dependence around the B 
mass region. 

First, if for two-body channels, 

M(s) = (1 - iksKk^M^s), (10) 

where M® are real, k is a diagonal momentum matrix and K a real-symmetric matrix, and 
the unitary S-matrix is given by [ |Tt| 

S = (l + ik^Kk^Xl-ik^Kk^)- 1 , (11) 

then it is easy to show that 

M = i(l + S)M (0) , (12) 
2 

for any number of two-body channels. 

Next we show that (12) ensures reality of the spectral function matrix <r(s) of (4). Using 
(10) and (12) in (4), we obtain 

a(s) = l(l-S+)(l + S)M°(s) 

= i(S-S+)M°(s), (13) 
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which is real. As S approaches a real-symmetric form, S + — > S, clearly the spectral function 
cr(s) approaches zero. Eq.(12) is not the only form that ensures reality of the spectral 
function. If we assume 

M(s) = S^M°(s), (14) 

one finds 

^) = ^-S^)M°(s) (15) 

which is also real. We comment later on (14). To see how (13) works and, in particular, 
to derive the results of f2|, let us specialize first to a 2-channel problem. The S-matrix is 
represented in terms of 3 real parameters (in general, an n-channel S matrix is expressed in 
terms of |(n + 1) real parameters), 

2^ ^2^gi(<5i+<52) ^g2i#2 

where 8\ and 5 2 are the eigenphases and rj, the elasticity, is limited by < rj < 1. Using 
(16) in (13), we obtain 

<tiM = \vsin25 lt2 Ml 2 {s) + ~(1 - V 2 )^cos(5 1 + 5 2 )M 2 ° 1 (s) (17) 

which are real. In addition, using (16) in (12), we get 

ImM lt2 {s) = ct 1i2 (s) (18) 
ReM h2 (s) = i(l + rjcos25 lj2 )Ml 2 - i(l - r ] 2 )hin(5 1 + 5 2 )M° ;1 . 

The result of that 

ImM^s) M$(s) 

ReM 1 (s) VC M 1 °( S ) 1 ' 

is recovered in the limit 5i j2 = and 77 = 1 — 2e with e small and positive. The argument of 
is that as inelasticity is needed to ensure reality of the spectral function, e can not vanish 
and Mi(s) is, therefore, complex implying that fsi effects do not disappear as the energy 
increases. In the following, we explore this issue further. 



This we do in the framework of a 3-channel problem, the reason being that apart from 
pedagogy, it allows for a more general form of the S-matrix — the off-diagonal elements of a 1- 
parameter representation of a 2-channel problem are purely imaginary while their 3-channel 
analogues are complex. 

The most general form of S-matrix in a 3-channel problem can be written as, 
Vie 2iSl (P12 + ^i 2 )e j(5l+<52) (pi 3 + ^is)e l{5l+<53) \ 

(Pi2 + iii2)e l ^ +5 ^ V 2e 2l&2 (p 23 + ^ 23 )e l(52+<53) (20) 

V (Pis + ^i3)e* (5l+ * 3) (p 23 + ^ 23 )e l(52+fe) r/ 3 e 2 ^ ) 

where with i ^ j ^ k, 

Pij = - -: * V ( 1 + + ^ ~ + ^ ~ ^ + Vj)( l -Vk + Vi- Vj)(l ~ Vk -Vi + Vj) 

&j = 7 ^_ V ( x + ^ ~ ^ ~ ^jX 1 + ^ + ^ + 77j)(l ~Vk + Vi + Vj)(Vk + Vi + Vj ~ 1) 

(21) 

The six real parameters are rji and Si(i = 1, 2, 3). In what follows we use (p + i£)ij = rjije 1 ^ . 
This unitary S-matrix has the following properties: First, if any one of rji becomes unity, 
the other two become identical. For example, if 773 = 1, then the requirement that all three 
off-diagonal 77^ be > implies that 771 = 772 = rj. In this situation channel 3 decouples 
from channels 1 and 2 and one recovers the 2-channel form (16) for the S-matrix for these 
channels. Second, if all 77^ — > 1, then all 77^ — > and the S-matrix becomes diagonal with 
diagonal entries e 2i<5 % i.e., each channel becomes elastic. Now using the S-matrix of (20) in 
(13) we can write down the analogues of (17) and (18) as follows: 

ff i( s ) = \ri x sin28 1 M^{s) + ^^sintyu + <&i + <5 2 )M 2 (s) 

+I m3S in(0 13 + S 1 + 5 3 )M 3 °(s) (22) 

ImMi(s) = <ti(s) 

ReM 1 (s) = ^(1 + r ]l cos25 1 )M°(s) + ^i 2 cos(0 12 + 5i + 5 2 )M° 2 {s) (23) 

+ ^iscos(0i3 + 5i + 5 3 )M 3 °(s). 



The analogous relations for other values of the index are written down by symmetry. Now, 
if the Pomeron were the sole contributor to the elastic amplitudes (the diagonal elements of 
the S-matrix), 5i would be strictly zero and rji close to ( and less than ) unity. The fact that 
Regge exchanges contribute to the elastic amplitudes makes Si finite through small . As 
\fs increases in the B mass region, the Regge-generated amplitudes move towards the base 
of the unitary circle while the Pomeron generated amplitudes do the same albeit at a much 
slower rate. Such a motion of the elastic amplitudes in the Argand diagram leads to an 
increase of rji towards unity and a decrease in the magnitudes of Si. 

A statement as to what happens as s — > oo depends very much on the behavior of 
the Pomeron amplitude. The Regge-exchange generated S-wave falls faster than s a (°) -1 . 
If, as s — > oo, the Pomeron "softens" such that ap(0) — > 1, then the Pomeron-generated 
S-wave amplitude will also move along the imaginary axis of the Argand diagram towards 
the base. Thus as s -> oo, all rji — > 1,%- — > and Si — * 0. The scattering amplitude 
reaches the elastic, though a "no-scattering" , limit. The S-matrix becomes a unit matrix and 
ImMi(s) -> 0, Mi(s) -> Mf If, on the other hand, as s — > oo the Pomeron-generated S- 
wave amplitudes reach fixed points on the imaginary axis, then all rji will approach values in 
the interval (0,1) with 2<5j = or 7r depending on whether the fixed point is below or above 
the center of the unitary circle. The S-matrix will approach a real-symmetric form. For 
three channels such an S-matrix has diagonal elements (t]i,t}2, tjs) and off-diagonal elements 
rjij with 

3 , 

Vi = 1 and T]^ = -^(1 - f]i)(l - vj). (24) 

i=i 

This form can generalized to an (nxn) case by constraining rji's by requiring Yh=i Vi = ( n ~ 2) 
and defining rjy as in (24). For the (3 x 3) case, if, say, rj 3 = 1, then rji = r]2 = and channel 
3 decouples from channels 1 and 2 which are now driven entirely by inelastic channels. If we 
assume that such a situation does not arise then rji < 1. This implies that elastic unitarity 
is not satisfied in any channel and inelastic channels must be there to implement unitarity. 
Eq.(13) now implies that <Ji{s) = (ImMi(s) = 0) while (12) implies that each of M;(s) 
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becomes a real superposition of all M-°\s). 

We conclude by a critique of the often-used adhoc procedure [fH[3yil| to unitarize real 
amplitudes M^(s) in presence of fsi given by (14). For a single channel, S = exp(2iS), and 
hence M(s) = exp(id)M°(s). This is said to be a statement of Watson's theorem. Eq.(14) 
is then a naive multi-channel generalization of the single- channel case. Watson's theorem 
, however, only claims that in a single-channel case, the phase of the formfactor, or decay 
amplitude, be the same as the scattering phase. It is to be emphasized that the magnitude 
of the formfactor also chang result of fsi. In general, for a single channel case the 



formfactors satisfy Muskhelishvili-Omnes equation fl2HI4 

Mis) = exp{- f 6 ^ dS ' }M°{s)exp{t5{s)) (25) 

71 J S — S 

where P signifies 'principal part'. Clearly, even for elastic scattering the magnitude of the 
formfactor changes - a circumstance assumption (14) forbids. 

In conclusion, we have investigated the issue of fsi in the limit of heavy quark masses 
using Pomeron and Regge exchanges as guides to high energy scattering. In the B mass 
region Pomeron exchange does not dominate as evidenced by the still-decreasing total cross- 
sections. As energy is increased ( but still in tens-of-GeV region), the S- matrix tends towards 
a real-symmetric form as the Regge-generated partial-wave amplitude drops with energy 
faster than s afi ' ' _1 . Asymptotically, if ap(0) remains at 1.08, partial-wave unitarity is 
violated around \/s ~ 10 u GeV. However, if the Pomeron were to "soften" such that 
op(0) — > 1 asymptotically then the S-matrix approaches a unit matrix-" no-scattering" 
limit. In this situation Mj(s) become real with Mj(s) — > M-°\s). If, on the other hand, the 
elastic amplitudes approach fixed points on the imaginary axis of the Argand diagram, the 
S-matrix becomes real-symmetric. Inelastic channels must be there to satisfy unitarity, yet 
they do not impart an imaginary part to the decay amplitude. Each M$(s) now becomes a 
real superposition of all M-;°'(s). 
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Sciences and Engineering Research Council of Canada to ANK. 
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